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We discuss an approach to obtaining black hole quasinormal modes (QNMs) using the asymptotic 
iteration method (AIM), initially developed to solve second order ordinary differential equations. 
We introduce the standard version of this method and present an improvement more suitable for 
numerical implementation. We demonstrate that the AIM can be used to find radial QNMs for 
Schwarzschild, Reissner-Nordstrom (RN) and Kerr black holes in a unified way. An advantage 
of the AIM over the standard continued fraction method (CFM) is that for differential equations 
with more than three regular singular points Gaussian eliminations are not required. However, the 
convergence of the AIM depends on the location of the radial or angular position, choosing the 
best such position in general remains an open problem. This review presents for the first time the 
spin 0, 1/2 & 2 QNMs of a Kerr black hole and the gravitational and electromagnetic QNMs of the 
RN black hole calculated via the AIM, and confirms results previously obtained using the CFM. 
We also presents some new results comparing the AIM to the WKB method. Finally we emphasize 
that the AIM is well suited to higher dimensional generalizations and we give an example of doubly 
rotating black holes. 
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I. INTRODUCTION 

The study of quasinormal modes (QNMs) of black holes is an old and well established subject, where the various 
frequencies are indicative of both the parameters of the black hole and the type of emissions possible. Initially the 
calculation of these frequencies was done in a purely numerical way, which requires selecting a value for the complex 
frequency, integrating the differential equation, and checking whether the boundary conditions are satisfied. Note 
that in the following we shall use the definition that QNMs are defined as solutions of the perturbed field equations 
with boundary conditions: 



. e x — > —oo 

a iuix X OO 



for an e~ lu>t time dependence (which corresponds to ingoing waves at the horizon and outgoing waves at infinity). 
Also note the boundary condition as x — > oo does not apply to asymptotically anti-de Sitter spacetimes, where 
instead something like a Dirichlet boundary condition is imposed, for example see Ref. [1]. Since those conditions 
are not satisfied in general, the complex frequency plane must be surveyed for discrete values that lead to QNMs. 
This technique is time consuming and cumbersome, making it difficult to systematically survey the QNMs for a wide 
range of parameter values. Following early work by Vishveshwara [2], Chandrasekhar and Detweiler [3] pioneered this 
method for studying QNMs. 

In order to improve on this, a few semi-analytic analyses were also attempted. In one approach, employed by 
Mashoon et al. [4], the potential barrier in the effective one-dimensional Schrodinger equation is replaced by a 
parameterized analytic potential barrier function for which simple exact solutions are known. The overall shape 
approximates that of the true black hole barrier, and the parameters of the barrier function are adjusted to fit the 
height and curvature of the true barrier at the peak. The resulting estimates for the QNM frequencies have been 
applied to the Schwarzschild, Reissner-Nordstrom and Kerr black holes, with agreement within a few percent with the 
numerical results of Chandrasekhar and Detweiler in the Schwarzschild case [3], and with Gunter [5] in the Reissner- 
Nordstrom case. However, as this method relies upon a specialized barrier function, there is no systematic way to 
estimate the errors or to improve the accuracy. 

The method by Leaver [6], which is a hybrid of the analytic and the numerical, successfully generates QNM 
frequencies by making use of an analytic infinite-series representation of the solutions, together with a numerical 
solution of an equation for the QNM frequencies which involves, typically by applying a Frobenius series solution 
approach, the use of continued fractions. This technique is known as the continued fraction method (CFM). 

Historically, another commonly applied technique is the WKB approximation [7]. Even though it is based on an 
approximation, this approach is powerful as the WKB approximation is known in many cases to be more accurate, 
and can be carried to higher orders, either as a means to improve accuracy or as a means to estimate the errors 
explicitly. Also it allows a more systematic study of QNMs than has been possible using outright numerical methods. 
The WKB approximation has since been extended to sixth-order [8]. 

However, all of these approaches have their limitations, where in recent years a new method has been developed 
which can be more efficient in some cases, called the asymptotic iteration method (AIM). Previously this method 
was used to solve eigenvalue problems [9] as a semi-analytic technique for solving second-order homogeneous linear 
differential equations. It has also been successfully shown by some of the current authors that the AIM is an efficient 
and accurate technique for calculating QNMs [10]. 

As such, we will review the AIM as applied to a variety of black hole spacetimes, making (where possible) compar- 
isons with the results calculated by the WKB method and the CFM a la Leaver [6]. Therefore, the structure of this 
paper shall be: In Sec. II we shall review the AIM and the improved method of Ciftci et al. [9] (also see Ref. [11]), 
along with a discussion of how the QNM boundary conditions are ensured. Applications to simple concrete examples, 
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such as the harmonic oscillator and the Poschl- Teller potential are also provided. In Sec. Ill the case of Schwarzschild 
(A)dS black holes shall be discussed, developing the integer and half-spin equations. In Sec. IV a review of how 
the QNMs of the Reissner-Nordstrom black holes shall be made, with several frequencies calculated in the AIM and 
compared with previous results. Sec. V will review the application of the AIM to Kerr black holes for spin 0, 1/2, 2 
fields. Sec. VI will discuss the spin-zero QNMs for doubly rotating black holes. We then summarize and conclude in 
Sec. VII. 

II. THE ASYMPTOTIC ITERATION METHOD 
A. The Method 

To begin we shall now review the idea behind the AIM, where we first consider the homogeneous linear second-order 
differential equation for the function xi x )i 

X " = X (x)x' + s (x) X , (2.1) 

where Ao(x) and Sq(x) are functions in C ao (a, b). In order to find a general solution to this equation, we rely on the 
symmetric structure of the right-hand of Eq. (2.1) [9]. If we differentiate Eq. (2.1) with respect to x, we find that 

X 1 " = Xi{x) X ' + si(x) X , 

where 

Ai = A„ + s + (Ao) 2 and si = s' + s Q X . 
Taking the second derivative of Eq. (2.1) we get 

X "" = A 2 (.*)*' + s 2 (x) X , 

where 

A2 = A' : + sx + A Ai and s x = s' + s A . 

Iteratively, for the (n + l) th and the (n + 2) th derivatives, n = 1, 2, we have 

X (n+1) = Xn-i(x) X ' + s n ^(x) X , (2.2) 

and thus bringing us to the crucial observation in the AIM is that differentiating the above equation n times with 
respect to x, leaves a symmetric form for the right hand side: 

X [n+2) = X n (x) X ' + s n (x)x , (2.3) 

where 

A„(ie) = A' n _ 1 (a;) + s n _ 1 (a:)+Ao(a;)A n _i(a;) and s n (x) = s'^^x) + s Q (x)X n _i(x) . (2.4) 

For sufficiently large n the asymptotic aspect of the "method" is introduced, that is: 

sj£) = s^(x) s 
X n {x) X n _i{x) 

where the QNMs are obtained from the "quantization condition" 

$n = s«A„-i - s n _iA„ = , (2.6) 
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which is equivalent to imposing a termination to the number of iterations [11]. From the ratio of the (n + \) th and 
the (n + 2) th derivatives, we have 



A v ("+2) 



A 



(x' + ftx) 

n-l (x' 



t^x 



(2.7) 



From our asymptotic limit, this reduces to 



dx 



m(x ( " +1) ) = 



An 

A n -i 



which yields 



X {n+1 \x) = d exp ^3)^') = CA,, 



_iexp / (^ + A )Gte' 



(2. 



(2.9) 



where C\ is the integration constant and the right-hand side of Eq. (2.4) and the definition of j3(x) have been used. 
Substituting this into Eq. (2.2), we obtain the first-order differential equation 



x' + /3x = Ciexp / ((3 + X )dx / 



which leads to the general solution 



X(x) = exp 



P(x')dx' 



C 2 + C! / exp 



[X (x") + 2/3(x")]dx"}dx' 



(2.10) 



(2.11) 



The integration constants, C\ and C2, can be determined by an appropriate choice of normalisation. Note, that for 
the generation of exact solutions C% = 0. 



B. The Improved Method 

Ciftci et al. [9] were among the first to note that an unappealing feature of the recursion relations in Eqs. (2.4) is 
that at each iteration one must take the derivative of the s and A terms of the previous iteration. This can slow the 
numerical implementation of the AIM down considerably and also lead to problems with numerical precision. 

To circumvent these issues we developed an improved version of the AIM which bypasses the need to take derivatives 
at each step [10]. This greatly improves both the accuracy and speed of the method. We expand the A n and s n in a 
Taylor series around the point at which the AIM is performed, £: 

00 

\ n (0 = J2<(x-0\ (2-12) 



=0 



where the c n and d l n are the i th Taylor coefficient's of A n (£) and s„(£) respectively. Substituting these expressions 
into Eqs. (2.4) leads to a set of recursion relations for the coefficients: 

< = (i + iRt 1 ! + dU + E c ° c ^ ' ( 2 - 14 ) 

i 

< = (i + lKt 1 + E4<-\- (2-15) 

k=0 
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In terms of these coefficients the "quantization condition" Eq. (2.6) can be re-expressed as 

- d<£ = , (2-16) 

and thus we have reduced the AIM into a set of recursion relations which no longer require derivative operators. 

Observing that the right hand side of Eqs. (2.14) and (2.15) involve terms of order at most n — 1, one can recurse 
these equations until only c and d l terms remain (that is, the coefficients of Ao and So only). However, for large 
numbers of iterations, due to the large number of terms, such expressions become impractical to compute. We avert 
this combinatorial problem by beginning at the n = stage and calculating the n + I coefficients sequentially until 
the desired number of recursions is reached. Since the quantisation condition only requires the i = term, at each 
iteration n we only need to determine coefficients with i < N — n, where N is the maximum number of iterations to 
be performed. The QNMs that we calculate in this paper will be determined using this improved AIM. 



2 + x 2 ) (f> = E<i). (2.17) 



C. Two Simple Examples 

1. The Harmonic Oscillator 

In order to understand the effectiveness of the AIM, it is appropriate to apply this method to a simple concrete 
problem: The harmonic oscillator potential in one dimension, 

£_ 

dx 2 

When \x\ approaches infinity, the wave function <j) must approach zero. Asymptotically the function <fr decays like a 
Gaussian distribution, in which case we can write 

0(z) = e~^f(x) , (2.18) 

where f(x) is the new wave function. Substituting Eq. (2.18) into Eq. (2.17) then re-arranging the equation and 
dividing by a common factor, one can obtain 

g=2*£ + (l- E)/ . (2.19) 

We recognise this as Hermite's equation. For convenience we let 1 — E = —2j, such that in our case Ao = 2x and 
so = — 2j. We define 

S n = A„s n _i - A n _is„ , for n= 1,2,3,... (2.20) 

Thus using Eqs. (2.4) one can find that 

n 

S n = 2" +1 [JO - i) , (2.21) 

i=0 

and the termination condition Eq. (2.6) can be written as 5 n = 0. Hence j must be a non-negative integer, which 
means 

E k = 2k + 1, for k = 0,1,2,... (2.22) 

and this is the exact spectrum for such a potential. Moreover, the wave function 4>(x) can also be derived in this 
method. 

We should point out that in this case the termination condition, 5 n = 0, is dependent only on the eigenvalue j for 
a given iteration number n, and this is the reason why we can obtain an exact eigenvalue. However, for the black 
hole cases in subsequent sections, the termination condition depends also on x, and therefore one can only obtain 
approximate eigenvalues by terminating the procedure after n iterations. 
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2. The Poschl-Teller Potential 



To conclude this section we will also demonstrate that the AIM can be applied to the case of QNMs, which have 
unbounded (scattering) like potentials, by recalling that we can find QNMs for Scarf II (upside-down Poschl- Teller- 
like) potentials [12]. This is based on observations made by one of the current authors [13] relating QNMs from 
quasi-exactly solvable models. Indeed bound state Poschl-Teller potentials have been used for QNM approximations 
previously by inverting black hole potentials [4]. However, the AIM does not require any inversion of the black hole 
potential as we shall show. 



Starting with the potential term 



and the Schrodinger equation, we obtain: 



dx 2 



V(x) = -sech 2 x , 



lo 2 — -sech x I tp = 



(2.23) 



(2.24) 



As we shall also see in the following sections, it is more convenient to transform our coordinates to a finite domain. 
Hence, we shall use the transformation y — tanhir, which leads to 



(i-y 2 ) 



dy 2 



dy 
■2y 



l-y 2 



dy 



(l-y 2 ) 2 2(1 -y 2 ) 



t/> = 0, 



(2.25) 



where — 1 < y < 1. The QNM boundary conditions in Eq. (1.1) can then be implemented as follows. As y — > 1 we 



shall have ip ~ e T«^ ^ (\ _ y)±W 2 . Hence our boundary condition ip ~ e 1 ' 



y — > — 1 we have if> 



(1 + y) ±lLU / 2 and the boundary condition ip ~ e 



can take the boundary conditions into account by writing 

^ = (l-y)-^(l + y)-^/ 2 4>, 



and therefore have 



where 



d 2 (j) _ 2y(l - iu>) d(j) 1 - 2iuj - 2uj 2 
dy 2 ~ l-y 2 dy + 2(1 - y 2 ) ' 



Ao 



2y(l - iu) 

l-y 2 ' 
1 - 2iu) - 2uj 2 



so = 



ip ~ (1 — y) lul / 2 . Likewise, as 
=4> V ~ (1 + y)~ iul/2 - As such we 

(2.26) 



(2.27) 

(2.28) 
(2.29) 



2(1 -V 2 ) ' 

Following the AIM procedure, that is, taking 8 n — successively for n = 1, 2, • • • , one can obtain exact eigenvalues: 



(2.30) 



This exact QNM spectrum is the same as the one in Ref. [13] obtained through algebraic means. 

The reader might wonder about approximate results for cases where Poschl-Teller approximations can be used, such 
as Schwarzschild and SdS backgrounds, e.g., see Refs. [1, 4]. In fact when the black hole potential can be modeled 
by a Scarf like potential the AIM can be used to find the eigenvalues exactly [12] and hence the QNMs numerically. 
We demonstrate this in the next section. 



III. SCHWARZSCHILD (A)DS BLACK HOLES 

We shall now begin the core focus of this review, the study of black hole QNMs using the AIM. Recall that the 
perturbations of the Schwarzschild black holes are described by the Regge- Wheeler [14] and Zerilli [15] equations, and 
the perturbations of Kerr black holes are described by the Teukolsky equations [16]. The perturbation equations for 
Reissner-Nordstrom black holes were also derived by Zerilli [17], and by Moncrief [18-20]. Their radial perturbation 
equations all have a one-dimensional Schrodinger-like form with an effective potential. 

Therefore, we shall commence in the coming subsections by describing the radial perturbation equations of 
Schwarzschild black holes first, where our perturbed metric shall be g^ = g ^ + h^, and where is spherically 
symmetric. As such it is natural to introduce a mode decomposition to h^. Typically we write 

* lm (t, r, e, 4>) = ^rue, 4>) , (3.1) 

r 

where Yi m (0,(j)) are the standard spherical harmonics. The function ui(r,t) then solves the wave equation 

^-u> 2 -Vt(rfj Ul (r)=0, (3.2) 

where x, defined by dx = dr/f(r), are the so-called tortoise coordinates and V(x) is a master potential of the form 
[21] 



V(r) = f(r) 

In this section 



(3.3) 



f(r) = l-^-K\ (3.4) 
r 3 



with cosmological constant A. Here s = 0,1,2 denotes the spin of the perturbation: scalar, electromagnetic and 
gravitational (for half-integer spin see Refs. [22-24] and Sec. VB). 

A. The Schwarschild Asymptotically Flat Case 

To explain the AIM we shall start with the simplest case of the radial component of a perturbation of the 
Schwarzschild metric outside the event horizon [15]. For an asymptotically flat Schwarzschild solution (A = 0) 

AO (3.5) 

where from dx = dr/f(r) we have 

x(r)=r + 2Mln(^-l) , (3.6) 

for the tortoise coordinate x. 

Note that for the Schwarzschild background the maximum of this potential, in terms of r, is given by [25] 



3M 1 



^+l)-(l-s 2 ) + (^(^+l) 2 + y^+l)(l- S 2 ) + (l- S 2 ) 2 ) 1/2 ] . (3.7) 



2 t(l+l) 

The choice of coordinates is somewhat arbitrary and in the next section (for SdS) we will see how an alternative 
choice leads to a simpler solution. Firstly, consider the change of variable: 

e = i-^, (3.8) 
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with < £ < 1. In terms of £, our radial equation then becomes 



d 2 ip l-3£ dtp 



+ 



4M 2 cj 2 £(£ + 1) 1 - s 2 



ip = 0. (3.9) 



To accommodate the out-going wave boundary condition ?/> — > e lu}x = e l "('"+ 2M in(r/2M-i)) ag ^ x ,r) — >• oo in terms 
of £ (which is the limit £ — > 1) and the regular singularity at the event horizon (£ — ¥ 0), we define 

m = r 2iM W(l _ f) -«M Ue ^ xK) ; (3 . 10 ) 

where the Coulomb power law is included in the asymptotic behaviour (cf. Ref. [6] Eq. (5)). The radial equation 
then takes the form: 

x" = M0x' + s (0x, (3-11) 

where 

, ( ,s _ 4A^(2g 2 -4£ + !)-(!- 3fl(l-fl 

A o(4) £(1 -C) 2 ' 

16Af 2 c 2 (g - 2) - 8Mia,(l - fl + ^ + !) + (!- * 2 )(1 - fl 

so(o = e(i-o 2 ' ( ^ 

Note that primes of x denote derivatives with respect to £. 

Using these expressions we have tabulated several QNM frequencies and compared them to the WKB method of 
Ref. [25] and the CFM of Ref. [6] in Table I. For completeness Table I also includes results from an approximate 
semi-analytic 3rd order WKB method [25]. More accurate semi-analytic results with better agreement to Leaver's 
method can be obtained by extending the WKB method to 6th order [8] and indeed in Sec. V we use this to compare 
with the AIM for results where the CFM has not been tabulated. 

It might also be worth mentioning that a different semi-analytic perturbative approach has recently been discussed 
by Dolan and Ottewill [26], which has the added benefit of easily being extended to any order in a perturbative 
scheme. 



The de-Sitter Case 



We have presented the QNMs for Schwarzchild gravitational perturbations in Table I, however, to further justify 
the use of this method, it is instructive to consider some more general cases. As such, we shall now consider the 
Schwarzschild de Sitter (SdS) case, where we have the same WKB-like wave equation and potential as in the radial 
equation earlier, though now 



2M , r 2 
/(r)=l- — -A T 



(3.14) 



where A > is the cosmological constant. Interestingly the choice of coordinates we use here leads to a simpler AIM 
solution, because there is no Coulomb power law tail; however, in the limit A = we recover the Schwarzschild results. 
Note that although it is possible to find an expression for the maximum of the potential in the radial equation, for 
the SdS case, it is the solution of a cubic equation, which for brevity we refrain from presenting here. In our AIM 
code we use a numerical routine to find the root to make the code more general. 

In the SdS case it is more convenient to change coordinates to £ = 1/r [1], which leads to the following master 
equation (cf. Eq. (3.3)) 



d 2 i> 
dt? 



p' dip 



£(£+1) + (1 - s 2 ) [2M£- (4-s 2 ) 



V = o, 



(3.15) 
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TABLE I. QNMs to 4 decimal places for gravitational perturbations (s — 2) where the fifth column is taken from Ref. [25]. 
Note that the imaginary part of the n = 0, I = 2 result in [25] has been corrected to agree with Ref. [6]. [*] Note also that 
if the number of iterations in the AIM is increased, to say 50, then we find agreement with Ref. [6] accurate to 6 significant 
figures. 



£ 


n 


^ Leaver 


maim (after 15 iterations) 


UWKB 


2 





0.3737 - 0.0896 i[*] 


0.3737 - 0.0896 i 


0.3732 - 0.0892 i 








(<0.01%)(<0.01%) 


(-0.13%) (0.44%) [*] 




1 


0.3467 - 0.2739 i 


0.3467 - 0.2739 i 


0.3460 - 0.2749 i 








(<0.01%)(<0.01%) 


(-0.20%) (-0.36%) 




2 


0.3011 - 0.4783 i 


0.3012 - 0.4785 i 


0.3029 - 0.4711 i 








(0.03%) (-0.04%) 


(0.60%)(1.5%) 




3 


0.2515 - 0.7051 i 


0.2523 - 0.7023 i 


0.2475 - 0.6703 i 








(0.32%) (0.40%) 


(-1.6%)(4.6%) 


3 





0.5994 - 0.0927 i 


0.5994 - 0.0927 i 


0.5993 - 0.0927 i 








(<0.01%)(<0.01%) 


(-0.02%) (0.0%) 




1 


0.5826 - 0.2813 i 


0.5826 - 0.2813 i 


0.5824 - 0.2814 i 








(<0.01%)(<0.01%) 


(-0.03%) (-0.04%) 




2 


0.5517 - 0.4791 i 


0.5517 - 0.4791 i 


0.5532 - 0.4767 i 








(<0.01%)(<0.01%) 


(0.27%) (0.50%) 




3 


0.5120 - 0.6903 i 


0.5120 - 0.6905 i 


0.5157 - 0.6774 i 








(<0.01%)(-0.03%) 


(0.72%)(1.9%) 




4 


0.4702 - 0.9156 i 


0.4715 - 0.9156 i 


0.4711 - 0.8815 i 








(0.28%)(<0.01%) 


(0.19%)(3.7%) 




5 


0.4314 - 1.152 i 


0.4360 - 1.147 i 


0.4189 - 1.088 i 








(1.07%)(0.43%) 


(-2.9%) (5.6%) 


4 





0.8092 - 0.0942 i 


0.8092 - 0.0942 i 


0.8091 - 0.0942 i 








(<0.01%)(<0.01%) 


(-0.01%)(0.0%) 




1 


0.7966 - 0.2843 i 


0.7966 - 0.2843 i 


0.7965 - 0.2844 i 








(<0.01%)(<0.01%) 


(-0.01%) (-0.04%) 




2 


0.7727 - 0.4799 i 


0.7727 - 0.4799 i 


0.7736 - 0.4790 i 








(<0.01%)(<0.01%) 


(0.12%)(0.19%) 




3 


0.7398 - 0.6839 i 


0.7398 - 0.6839 i 


0.7433 - 0.6783 i 








(<0.01%)(<0.01%) 


(0.47%) (0.82%) 




4 


0.7015 - 0.8982 i 


0.7014 - 0.8985 i 


0.7072 - 0.8813 i 








(-0.01%) (-0.03%) 


(0.81%)(1.9%) 



where we have defined 

V = i 2 - 2M£ 3 - A/3 =>■ p' = 2£(1 - 3M£) . (3.16) 

It may be worth mentioning that for SdS we can express [1] : 

e iux = (£ - fr)^T(£ - - (3.17) 

in terms of the roots of f(r), where £i is the event horizon and £2 is the cosmological horizon (and K n is the surface 
gravity at each £„). This is useful for choosing the appropriate scaling behaviour for QNM boundary conditions. 
Based on the above equation an appropriate choice for QNMs is to scale out the divergent behaviour at the 
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cosmological horizon: 1 



which implies 



pu" + {p' - 2iui)u' 



t(l + 1) + (1 - s^) 2M£ - (4 - s z ) 



A 
'6^ 



71 = , 



(3.18) 



(3.19) 



in terms of £. Furthermore, based on the scaling in Eq. (3.18), the correct QNM condition at the horizon £i implies 

= %XW, (3.20) 

where 



1# 

2dr 



= _ 



1 A 
36 



(3.21) 



with £i — 1/ri, and n is the smallest real solution of /(r) = 0, implying p = 0. The differential equation then takes 
the standard AIM form: 



x" = Ao(Ox' + *o(Ox, 



where 



Ao(0 
*o(0 



2iuj 



2iu> 



£(£- 



l) + (l-s 2 ) ( 2M£- (4-s 2 )- 



A 



/iuj 



1 



(p' - 2iw) 



(3.22) 

(3.23) 
(3.24) 



'6£V ^i(£-£i) 2 W 
Using these equations, we present in Table II results for SdS with s = 2. 

Identical results were generated by the AIM and CFM, both after 50 iterations. Though results are presented for 
n = 1, 2, 3, t = 2, 3 modes only, the AIM is robust enough to be applied to any other case; where like the A = case, 
agreement with other methods in more extreme parameter choices would only require further iterations. 

As far as we are aware only Ref. [27] (who used a semi-analytic WKB approach) has presented tables for general 
spin fields for the SdS case. We have also compared our results to those in Ref. [27] for the s = 0, 1 cases and find 
identical results (to a given accuracy in the WKB method) . 

It may be worth mentioning that a set of three-term recurrence relations was derived in Ref. [10] for the CFM, 
valid for electromagnetic and gravitational perturbations (s = 1,2), while for s = this reduces to a five-term 
recurrence relation. However, for the AIM we can treat the s — 0,1,2 perturbations on an equal footing, see Ref. 
[10] for more details. Typically n Gaussian Elimination steps are required to reduce an n + 3 recurrence to a 3-term 
continued fraction, e.g., for Reissner-Nordtrom see Ref. [28] and for higher-dimensional Schwarzschild backgrounds 
see Ref. [29] (for an application of the CFM to higher dimensional asymptotic QNMs see [30]). However, all that is 
necessary in the AIM is to factor out the correct asymptotic behaviour at the horizon(s) and infinity (we showed this 
for higher-dimensional scalar spheroids in Ref. [31]). 



C. The Spin-Zero Anti-de-Sitter Case 

There are various approaches to finding QNMs for the SAdS case (an eloquent discussion is given in the appendix 
of Ref. [32], see also [33]). One approach is that of Horowitz and Hubeny [34], which uses a series solution chosen 



1 Note that this is opposite to the case presented in Ref. [1], where they define the QNMs as solutions with boundary conditions 
ip(x) oc e^ LlLJX as x — > ±oo, for e 1 "* time dependence. 
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TABLE II. QNMs to 6 significant figures for Schwarzschild de Sitter gravitational perturbations (s — 2) for 1 = 2 and 1 = 3 
modes. We only present results for the AIM method, because the results are identical to those of the CFM after a given number 
of iterations (in this case 50 iterations for both methods). The n — 1,2 modes can be compared with the results in Ref. [27] 
for s = 2. 



A(£ = 


2) 


n 


= 1 


n 


= 2 


n = 


= 3 







0.373672 - 


0.0889623 i 


0.346711 


- 0.273915 i 


0.301050 - 


0.478281 i 






U.OOOOiJl- - 


fl 081 7^84 \ 


fl 31 87^0 
u.oio i oy 


n 94Q1 07 i 


n 989739 


n 49Q484 i 






n 90880^ 


07390R7 i 
u.u i ozyu i i 


fl 98^841 


n 991 794 i 


n 9^0009 
u.zuyyyz - 


n 377009 i 

u.o i i uyz i 


OR 
u.uu 




9 c i' : i98Q 


n 08^49^ i 


fl 94^749 


n 1 8Q7Q1 i 
u.ioy i yi i 


n 93007fi 


31 Q1 ^7 i 


n ns 

u.uo 




n 1 Q7489 


fl n40877' I i i 


1 Q41 1 ^ 

u. ±y^±± lo 


n 1 4Q787 i 
- u. ±^±y jo I i 


n 1 871 9n 

U.IOI ±zu - 


n 9^09^7 i 


n no 




n i R9Ri n 


U.UIlOUuU 1 


1 R0780 
u. iuu i oy 


n 1 941 ^9 i 


n i c i7n49 

U. lO 1 U1Z - 


n 9n7i 1 7 i 

U.iUI 111 1 


n in 

U.1U 




nil 7Q1 fi 


n n^n9i n^ i 


fill 7943 


n oonfi4no i 


nil ^R7R 


n i ^i i n9 i 


n 1 1 




u.uo i iuyy 


n nnoisi ^fi^ i 


0379403 
u.uo i i'lyo 


n n98847n i 


0379081 

U.UO 1 iUOl 


n n480784 i 


A {£ = 


3) 


n 


= 1 


n 


= 2 


n - 


= 3 







0.599443 - 


0.0927030 i 


0.582644 


- 0.281298 i 


0.551685 - 


0.479093 i 


0.02 




0.543115 - 


0.0844957 i 


0.530744 


- 0.255363 i 


0.507015 - 


0.432059 i 


0.04 




0.480058 - 


0.0751464 i 


0.471658 


- 0.226395 i 


0.455011 - 


0.380773 i 


0.06 




0.407175 - 


0.0641396 i 


0.402171 


- 0.192807 i 


0.392053 - 


0.322769 i 


0.08 




0.317805 - 


0.0503821 i 


0.315495 


- 0.151249 i 


0.310803 - 


0.252450 i 


0.09 




0.261843 - 


0.0416439 i 


0.260572 


- 0.124969 i 


0.257998 - 


0.208412 i 


0.10 




0.189994 - 


0.0303145 i 


0.189517- 


0.0909507 i 


0.188555 - 


0.151609 i 


0.11 




0.0600915 - 


0.00961888 i 


0.0600766 


- 0.0288567 i 


0.0600469 - 


0.0480945 i 



to satisfy the SAdS QNM boundary conditions. This method can easily be applied to all perturbations (s = 0, 1, 2). 
The other approach is to use the Frobenius method of Leaver [6] , but instead of developing a continued fraction the 
series must satisfy a boundary condition at infinity, such as a Dirichlet boundary condition [1]. 

The AIM does not seem easy to apply to metrics where there is an asymptotically anti-de Sitter background, 
because for general spin, s, the potential at infinity is a constant and hence would include a combination of ingoing 
and outgoing waves, leading to a sinusoidal dependence [35]. However, for the scalar spin zero (s = 0) case, the 
potential actually blows up at infinity and is effectively a bound state problem. In this case the AIM can easily be 
applied as we show below. 

Let us consider the scalar wave equation in SAdS spacetime, where A = — 3/i? 2 , and R is the AdS radius. The 
master equation takes the same form as for the graviational case, except that the potential becomes 

v ={ i ~l+ '-') + 2 ) = 2(r ~ ~ V* + 2><r ' + 11 • ^ 

Here for simplicity we have taken the AdS radius R = 1, the mass of the black hole M = 1, and the angular momentum 
number I — 0. Hence the horizon radius r + = 1. Thus, with this choice we can compare with the data in Table 3.2 
on page 37 of Ref. [36] (see Table III below). 

To implement the AIM we first look at the asymptotic behavior of ip. As r — s- r + = 1, the potential V goes to zero. 
In addition, 
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For QNMs we choose the "out-going" (into the black hole) boundary condition. That is, 



$ ~ e~ wx ~ I 1 - - 1 . (3.27) 

On the other extreme of our space, r — ¥ oo, the potential goes to infinity. This is a crucial difference from the case 
of gravitational perturbations. In that case, the potential goes to a constant. However, in the scalar case, as r — > oo, 
if) ~ (l/r) ±v/2+1 ^ 2 and to implement the Dirichlet boundary condition, we take 

V>~(-J . (3.28) 



For the AIM one possible choice of variables is 



£ = 1 - 1 , (3.29) 
r 



and we see that to accommodate the asymptotic behaviour of the wavefunction we should take 

V-> = r W4 (l-0^X- (3-30) 
Finally, after some work we find the scalar perturbation equation is 

x" = M0x' + s (0x, (3-31) 

where 

-iuq + 2[-4 + 2(9 + 4V2)£, - (21 + 10^2)^ + 4(2 + V2)£ 3 } 



An = -- 



2£q 



s = — ^— Uiu[9 + 8V2- 2(7 + 5^)£ + (6 + <±V2)£ 2 ]q + w 2 (-l + £) 2 (-40 + 41£ - 20£ 2 + 4£ 3 
-4[4 - 5£ + 2^ 2 ][-8(3 + 2V2) + 8(10 + 7^)£ - (91 + 64\/2)£ 2 + (34 + 24V2)£ 3 ] J , 

and q = (-4 + 9£ - 7£ 2 + 2£ 3 ). Using the AIM we find the results presented in Table III below. 



(3.32) 



(3.33) 



TABLE III. Comparison of the Srst few QNMs to 6 significant figures for Schwarzschild anti de Sitter scalar perturbations 
(s — 0) for £ = modes with r+ = 1. The second column corresponds to data [36] using the Horowitz and Hubeny (HH) 
method [34], while the third column is for the AIM using 70 iterations. [*] Note the mismatch for the real part of the n = 3 
mode in Ref. [36]; we have confirmed this using the Mathematica notebook provided in Ref. [21]. 



n 


HH method 


AIM 





2.7982 - 2.6712 i 


2.79823 -2.67121 i 


1 


4.75849 - 5.03757 i 


4.75850 -5.03757 i 


2 


6.71927 - 7.39449 i 


6.71931 -7.39450 i 


3 


8.68223[*] - 9.74852 i 


8.68233 -9.74854 i 


4 


10.6467 - 12.1012 i 


10.6469 -12.1013 i 


5 


12.6121 - 14.4533 i 


12.6125 -14.4533 i 


6 


14.5782 - 16.8049 i 


14.5788 -16.8050 i 


7 


16.5449 - 19.1562 i 


16.5457 -19.1563 i 
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IV. REISSNER-NORDSTROM BLACK HOLES 



The procedure for obtaining the quasinormal frequencies of Reissner-Nordstrom black holes in four-dimensional 
spacetime is similar to that of our earlier cases. Starting with the Reissner-Nordstrom metric 



ds 2 



-f{r)dt 2 



dr 2 
W) 



r 2 d9 2 



r 2 sin 2 



(4.1) 



where f(r) = \1 — i + and \Q\ < | is the charge of the black hole. If we consider perturbations exterior to 
the event horizon, the perturbation equations of the Reissner-Nordstrom (charged and non-rotating) geometry can be 
separated into two pairs of Schrodinger-like equations, which describe the even- and odd-parity oscillations respectively 
[17-20]. They are given by 



dx 2 



o 



where (+) corresponds to even- and (- 



Vt\r) 



to odd-parity modes: 
4Q 2 



A 



Ar - 



9j 



V^\r) = vt\r) + 2 qj d 



dx \r 2 [(l-l){l + 2)r + q/ 



for i = j = 1,2 (i =/= j) where 



dr 
dx 
A 

A 



A 

r 

2 



1(1- 

1 r 



+ Q 2 = (r-r+)(r-r-) 



qi = 2 



?2 



1 r 

2 



v/9 + 16Q 2 (Z- l)(Z + 2) 
v/9 + 16Q 2 (Z- l)(Z + 2) 



(4.2) 

(4.3) 
(4.4) 

(4.5) 
(4.6) 
(4.7) 
(4.8) 

(4.9) 



and p — —ioj. Here u is the frequency, I the angular momentum parameter and r_ and r + the radii of the inner and 
outer (event) horizons of the black hole respectively. Note that r + — 1 and r_ = at the Schwarzschild limit (Q — 0); 



2 at the extremal limit (Q — \). Here the tortoise coordinate is given by 



-dr 



ln(r — r_|_) 



ln(r 



(4.10) 



which ranges from — oo at the event horizon to +oo at spatial infinity. 

The QNMs of the Reissner-Nordstrom black holes are ordinarily accompanied by the emission of both electro- 
magnetic and gravitational radiation, except at the Schwarzschild limit [28, 37]. Eq. (4.2) corresponds to purely 
gravitational perturbations for the radial wave functions Z% and purely electromagnetic perturbations for Z± at 



the Schwarzschild limit. Chandrasekhar 

A-) 



has shown that the solution for the even-parity oscillations and 



Z\ for the odd-parity oscillations have the relationship 



[A(A - 2) - 2pq j ] Z 



( + ) 



A(A-2) 



2q 2 A 



r 3 [(A - 2)r + q,] 



Z 



(-) 



2qy 



dZ 



(-) 



dx 



(4.11) 



so one can just consider solutions for a specific parity, as in the Schwarzschild case, to understand the property of the 
black hole. Since the formalism of the effective potential ^ in the odd-parity equation is much simpler than 
in the even-parity equation, it is customary to compute the QNMs for the odd-parity modes. 
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Note that the mass M of the Reissner-Nordstrom black hole has been scaled to 2M = 1, so its quasinormal 
frequencies are uniquely determined by the charge Q, the angular momentum I, and the overtone number n of the 
mode. 

The following procedure is similar to that in Sec. III. At first we change r to the variable x in Eq. (4.2) for the 
odd-parity mode. From Eq. (4.5) we have 



and 



dx 2 



d _ A d 
dx r 2 dr ' 



A\ ft - 2Q 2 \ d 
r 2 1 V r 3 ) dr 



A\ 2 d?_ 
dr 2 



Substituting Eq. (4.13) into Eq. (4.2) for the odd-parity mode, we get 



2 d^ 
dr 2 



Z 



A\ fr-2Q 2 



dr 



(-) 



P 2 + V i 



(-) 



z 



(-) 



Considering the QNM boundary conditions in the Reissner-Nordstrom case 



. 



(4.12) 



(4.13) 



(4.14) 



z {-) I e lulx ;x^-oo 



(4.15) 



and incorporating this into the radial wave function z\ ' , we have a form involving the asymptotic behaviour [28] 



Z 



(-) 



pr r 1 (r — r-) 1 P r + r - (r — r + ) r + T - XZi(r) 



(4.16) 



Differentiating Eq. (4.16) one and two times with respect to r, we have 



Z 



(-) 



,(-) 



-pr„,-l 



dr 



and 



(4.17) 



'(-) _ B -pr„-l t 



(r-r-) 1 9 (r-r + y+-'- [ XZi ,rr + ^zXz^r + (r| + T Z ,r) Xz % ] 



where Tz is defined by 



T z = -p- 



1,(1 P)(r+ - r-) - pr\ 



pr\ 



(r_|_ — r_)(r — r_) (r + — r_)(r — r + ) 
Substituting Eqs. (4.17) and (4.18) into Eq. (4.14), we obtain 

'A\ 2 /A\ /r-2Q 2 ' " 



A s 

Z2 ) XZt.n- 



2T t 



Xz it r 



(r| + r z , r ) + 



A\ fr-2Q 2 



p 2 + vt } 



Xz % =0 



(4.18) 



(4.19) 



(4.20) 



For the same reason as in Sec. Ill, here we change the variable r to £ by the definition £ = 1 — ^r, which ranges 
from at the event horizon to 1 at spatial infinity. Thus we have 



d 
dr 



(W) 2 d 



(4.21) 
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and 



d^ 
dr 2 



(i-0 4 d 2 _ 2 (lr0 3 d 



r\ d£ 2 r% di 

Substituting Eqs. (4.21) and (4.22) into Eq. (4.20), and rewriting the equation in the AIM form, we obtain 



where 



szM) 



2 



2r+T z 

a-e) 2 



-2Q 2 (1-Q 
A(l-€) 2 



Tz = -p- 



A»(l-0 8 A(1-£) 4L 
1-e ( [(l-p)(r+-r_ 
r+ (r+-r_)[r + 



2Q 2 (l-0] 

21 (i-o 



dr 2 



(-) 



= A 



A = 



(1 - £) 5 [ Ar+ 
[r + -r_(l-0] 



4Q 2 (1 



— pr 
r_(l 

-0 



0] 



A(l-0 4 

pr+(l-Q 

(r+-r_)€ 



(i-0 s 



(4.22) 

(4.23) 

(4.24) 

(4.25) 
(4.26) 
(4.27) 
(4.28) 



(i-0 2 

The numerical results to four decimal places are presented in Table IV, V and VI. They are compared with p Leaver 
and pwkb from Refs. [28] and [37] respectively. The quasinormal frequencies appear as complex conjugate pairs in 
p; we list only the ones with Im(p) > 0. Note that we arrange p as (Im(p) , Re(p)) . In Table VI the quasinormal 
frequencies obtained by the WKB method are not available. It is apparent that the quasinormal frequencies obtained 
by the AIM are very accurate except for n = 2 in the extremal case Q = ^ in Tables IV and VI. 

The QNMs of I = 2 and i = 2 in Table IV reduce to the purely gravitational QNMs in the Schwarzschild case at 
Q = 0, while the QNMs of I — 2 and i = 1 in Table V reduce to the purely electromagnetic QNMs at Q = 0. 

Some comments on the higher (n = 2) overtones for the Reissner-Nordstrom black hole for the extremal limit 
(Q = 1/2) are perhaps necessary. In general, much like the CFM the AIM begins to break down for larger overtones, 
requiring more iterations. However, near the extremal limit (Q = 1/2) the horizons become degenerate and the 
singularity structure of the corresponding differential (radial) equation changes [39] (the number of singular points 
are different in the non-extremal and the extremal cases), and causes the current implementation of the AIM (cf. 
4.16) to break down. Thus, we see in Tables IV- VI that for Q = 0.495 some of the values have large errors when 
compared to the CFM. 



V. KERR BLACK HOLES 



A rotating black hole carrying angular momentum is described by the Kerr metric (in Boyer-Lindquist coordinates) 



as 



with 



ds 2 



1 - - ] 



2ar sin 2 8 



dtdcj) + —dr 2 Y;d6 2 + [r 2 + a 2 



A = r 2 + a 2 - 2Mr = (r - r_)(r - r+) 
S = r 2 + a 2 cos 2 9 , 



a 2 r sin 2 9 



(5.1) 



(5.2) 
(5.3) 



and where a is the Kerr rotation parameter with < a < M, M being included as a general black hole mass. The 
horizons r_ and r + are again the inner and the outer (event) horizons respectively. Teukolsky [16] showed that 
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TABLE IV. Reissner-Nordstrom quasinormal frequency parameter values (p = — icu) for the fundamental (n — 0) and two 
lowest overtones for I = 2 and i — 2. 



n 


Q 


pLeaver 


pAIM 


PWKB 








(0.7473,-0.1779) 


(0.7473,-0.1779) 


(0.7463,-0.1784) 








(<0.01%)(<0.01%) 


(-0.13%) (-0.28%) 





0.2 


(0.7569,-0.1788) 


(0.7569,-0.1788) 


(0.7558,-0.1793) 








(<0.01%)(<0.01%) 


(-0.15%) (-0.28%) 





0.4 


(0.8024,-0.1793) 


(0.8024,-0.1793) 


(0.8011,-0.1797) 








(<0.01%)(<0.01%) 


(-0.16%) (-0.22%) 





0.495 


(0.8586,-0.1685) 


(0.8586,-0.1685) 


(0.8566,-0.1706) 








(<0.01%)(<0.01%) 


(-0.23%) (-1.25%) 


1 





(0.6934,-0.5478) 


(0.6934,-0.5478) 


(0.6920,-0.5478) 








(<0.01%)(<0.01%) 


(-0.20%) (-0.37%) 


1 


0.2 


(0.7035,-0.5503) 


(0.7035,-0.5502) 


(0.7020,-0.5522) 








(<0.01%)(0.02%) 


(-0.21%) (-0.36%) 


1 


0.4 


(0.7538,-0.5499) 


(0.7538,-0.5499) 


(0.7510,-0.5525) 








(<0.01%)(<0.01%) 


(-0.37%) (-0.47%) 


1 


0.495 


(0.8070,-0.5140) 


(0.8067,-0.5164) 


(0.8068,-0.5287) 








(-0.04%) (0.47%) 


(0.01%) (-2.21%) 


2 





(0.6021,-0.9566) 


(0.6021,-0.9566) 


(0.6059,-0.9421) 








(<0.01%)(<0.01%) 


(0.63%) (1.52%) 


2 


0.2 


(0.6129,-0.9599) 


(0.6128,-0.9599) 


(0.6164,-0.9458) 








(0.02%)(<0.01%) 


(0.57%) (1.47%) 


2 


0.4 


(0.6703,-0.9531) 


(0.6703,-0.9531) 


(0.6717,-0.9455) 








(<0.01%)(<0.01%) 


(0.21%) (0.80%) 


2 


0.495 


(0.7078,-0.8872) 


(0.8350,-0.8347) 


(0.7344,-0.9135) 








(17.97%) (5.92%) 


(2.66%) (-2.96%) 



the perturbation equations in the Kerr geometry are separable, where the separated equations for the angular wave 
function s Si m (9) and the radial wave function R(r) are given by: 



{(l-u 2 )S, u ], u + 



a 2 ui 2 u 2 — 2aujsu + s + s Ai m 



(m + su) z 



l~u 2 

AR. rr + (s + l)(2r - l)R r + K(r)R = 



(5.4) 
(5.5) 



where the function 



K(r) = — | (r 2 + a 2 ) 2 u? - 2amojr + a 2 m 2 + is [am (2r - 1) - u> (r 2 - a 2 )] j + 2is 



lm 



(5.6) 



In the above u = cos#, s is the spin weight, s Ai m is the spin-weighted separation constant for the angular equation, 
and m is another angular momentum parameter. For completeness the evaluation of the separation constant s Ai m 
using the AIM is discussed in Appendix A. 

In order to use the AIM we need to solve for the angular solution in the radial equation. However, for nonzero s 
the effective potential of the radial equation is in general complex. A straight forward application of the AIM does 
not give the correct answer. In fact a similar problem occurs in both numerical [40] and WKB [7] methods. For this 
reason we shall look at each of the spin cases (0, ^, 2) separately in the following subsections. 
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TABLE V. Reissner-Nordstrom quasinormal frequency parameter values (p — —ito) for the fundamental (n — 1) and two lowest 
overtones for I = 2 and i = 1. 



n 


Q 


pLeaver 


pAIM 


PWKB 








(0.9152,-0.1900) 


(0.9152,-0.1900) 


(0.9143,-0.1901) 








(<0.01%)(<0.01%) 


(-0.10%) (-0.05%) 





0.2 


(0.9599,-0.1929) 


(0.9599,-0.1929) 


(0.9590,-0.1930) 








(<0.01%)(<0.01%) 


(-0.09%) (-0.05%) 





0.4 


(1.1403,-0.1984) 


(1.1403,-0.1981) 


(1.1395,-0.1980) 








(<0.01%)(0.15%) 


(-0.07%) (-0.20%) 





0.495 


(1.3855,-0.1773) 


(1.3855,-0.1773) 


(1.3850,-0.1783) 








(<0.01%)(<0.01%) 


(-0.04%) (-0.56%) 


1 





(0.8731,-0.5814) 


(0.8731,-0.5814) 


(0.8717,-0.5819) 








(<0.01%)(<0.01%) 


(-0.16%) (-0.09%) 


1 


0.2 


(0.9200,-0.5894) 


(0.9200,-0.5894) 


(0.9186,-0.5897) 








(<0.01%)(<0.01%) 


(-0.15%) (-0.05%) 


1 


0.4 


(1.1100,-0.6021) 


(1.1100,-0.6021) 


(1.1081,-0.6014) 








(<0.01%)(<0.01%) 


(-0.17%)(0.12%) 


1 


0.495 


(1.3573,-0.5350) 


(1.3573,-0.5350) 


(1.3579,-0.5423) 








(<0.01%)(<0.01%) 


(0.04%) (-1.36%) 


2 





(0.8024,-1.0032) 


(0.8024,-1.0032) 


(0.8046,-0.9917) 








(<0.01%)(<0.01%) 


(0.27%) (1.15%) 


2 


0.2 


(0.8530,-1.0143) 


(0.8530,-1.0143) 


(0.8548,-1.0037) 








(<0.01%)(<0.01%) 


(0.21%) (1.05%) 


2 


0.4 


(1.0582,-1.0263) 


(1.0582,-1.0263) 


(1.0568,-1.0181) 








(<0.01%)(<0.01%) 


(-0.13%) (0.80%) 


2 


0.495 


(1.3019,-0.9024) 


(1.3019,-0.9024) 


(1.3141,-0.9222) 








(<0.01%)(<0.01%) 


(0.94%) (-2.19%) 



A. The Spin-Zero Case 

Because the AIM works better on a compact domain, we define a new variable y = 1 — which ranges from at 
the event horizon (r = r_|_) to 1 at spatial infinity. It is then necessary to incorporate the boundary conditions, which 
expressed in the new compact domain, where 



R(r) = (r 2 + a 2 )- 1 / 2 ^) 



(5.7) 



is 

By making the change of coordinates and change of function, Eq. (5.8) takes the form 

x(y) = A (y) + so(y) , 

where 

„ldff Id/ 



An 



'gdy f dy ' 



(5.8) 



(5.9) 



(5.10) 
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TABLE VI. Reissner-Nordstrom quasinormal frequency parameter values (p = — icu) for the fundamental (n — 1) and two 
lowest overtones for I = 1 and i — 1. 



n 


Q 


pLeaver 


pAIM 








(0.4965,-0.1850) 


(0.4965,-0.1850) 








(<0.01%)(<0.01%) 





0.2 


(0.5238,-0.1883) 


(0.5238,-0.1883) 








(<0.01%)(<0.01%) 





0.4 


(0.6470,-0.1965) 


(0.6470,-0.1965) 








(<0.01%)(<0.01%) 





0.495 


(0.8428,-0.1742) 


(0.8428,-0.1742) 








(<0.01%)(<0.01%) 


1 





(0.4290,-0.5873) 


(0.4292,-0.5873) 








(<0.01%)(0.05%) 


1 


0.2 


(0.4598,-0.5953) 


(0.4598,-0.5953) 








(<0.01%)(<0.01%) 


1 


0.4 


(0.5980,-0.6107) 


(0.5980,-0.6107) 








(<0.01%)(<0.01%) 


1 


0.495 


(0.7979,-0.5293) 


(0.7978,-0.5280) 








(0.01%) (0.25%) 


2 





(0.3496,-1.0504) 


(0.3495,-1.0504) 








(0.03%)(<0.01%) 


2 


0.2 


(0.3832,-1.0596) 


(0.3832,-1.0596) 








(<0.01%)(<0.01%) 


2 


0.4 


(0.5340,-1.0660) 


(0.5340,-1.0659) 








(<0.01%)(<0.01%) 


2 


0.495 


(0.7104,-0.9055) 


(0.6248,-1.0574) 








(-12.05%) (-16.78%) 



and 



1 d 2 .g 1 df 1 dg 1 , 2 , 



In the above we have defined 



and 



A dy 

r 2 + a 2 dr 



(5-11) 



(5.12) 



-»o-+ >r + (l-j) 1 



(5.13) 



where A = r l + a 2 — 2Mr, 



o-± 



-[(r± + a 2 )uj + am] , 



(5.14) 



and a is again our rotation parameter. Eq. (5.9) is now in the correct form to use the AIM for QNM frequency 
calculations. Note that the potential, V is 

1 



V = : --((X-2(r 2 + fl 2 )) 2 -AA) 



(5.15) 
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where the angular separation constant is defined via A = A^ m + a 2 u! 2 — 2amuj. 2 

As presented in Tables VII and VIII are the QNM frequencies for the scalar perturbations of the Kerr black hole 
with the two "extreme" (minimum and maximum) values of the angular momentum per unit mass, that is, a = 0.00 
and a = 0.80. m was set to 0, while I was given values of 0, 1 and 2 and n varied accordingly. 

Included in Table VII are the numerically determined QNM frequencies published by Leaver in 1985 [6]. The 
percentages bracketed under each QNM frequency via the AIM are the percentage differences between the calculated 
value and the numerical value published by Leaver. With the exceptions of the QNM frequencies for I = 0, n = and 
I = 2, n = 2, the AIM values correspond to the CFM up to four decimal places and even those "anomalies" differ by 
less than 0.30%. Proving, at least in this case, the AIM is a precise semi-analytical technique. 

In Table VIII, all three values were calculated in this work, even though published values are available for the third 
order WKB(J), at least graphically, where numerical values using the CFM were taken from Ref. [41]. Since the 
WKB(J) is a generally accepted semi- analytical technique for QNM frequency calculations, the percentages below the 
AIM values are the differences to the sixth order WKB(J) values. Only in the case of I — 0, n = does the AIM QNM 
frequency significantly differ from the sixth order WKB(J) value. Note that in an upcoming work, further values will 
be presented for values of a = 0.20, a = 0.40 and a = 0.60, with the same variations of I and n with M = 1 [42]. 



TABLE VII. The QNM Frequencies for the Scalar Perturbations of the Kerr Black Hole, with a — 0.00, that is, the Schwarzschild 
limit (M = 1, m — 0). Numerical data via the CFM taken from [41], where the AIM was set to run at 15 iterations. 



1 


n 


Numerical 


Third Order WKB(J) 


Sixth Order WKB(J) 


AIM 








0.1105 - 0.1049i 


0.1046 - O.lOOli 


0.1105 - 0.1008i 


0.1103 - 0.1046i 








(-5.34% , 9.82%) 


(<0.01% , -2.91%) 


(-0.18% , -0.29%) 


1 





0.2929 - 0.0977i 


0.2911 - 0.0989i 


0.2929 - 0.0978i 


0.2929 - 0.0977i 








(-0.61% , 1.23%) 


(<0.01% , 0.10%) 


(<0.01% , <0.01%) 




1 


0.2645 - 0.3063i 


0.2622 - 0.3074i 


0.2645 - 0.3065i 


0.2645 - 0.3063i 








(-0.87% , 0.36%) 


(<0.01% , 0.07%) 


(<0.01% , <0.01%) 


2 





0.4836 - 0.0968i 


0.4832 - 0.0968i 


0.4836 - 0.0968i 


0.4836 - 0.0968i 








(-0.08% , <0.01%) 


(<0.01% , <0.01%) 


(<0.01% , <0.01%) 




1 


0.4639 - 0.2956i 


0.4632 - 0.2958i 


0.4638 - 0.2956i 


0.4639 - 0.2956i 








(-0.15% , 0.07%) 


(-0.02% , <0.01%) 


(<0.01% , <0.01%) 




2 


0.4305 - 0.5086i 


0.4317 - 0.5034i 


0.4304 - 0.5087i 


0.4306 - 0.5086i 








(0.28% , -1.02%) 


(-0.02% , 0.02%) 


(0.02% , <0.01%) 



B. The Spin-Half Case 

For the spin-1/2 case we would like to know how the AIM can be used to derive an appropriate form of the Dirac 
equation in this spacetime background using the basis set up by four null vectors which are the basis of the Newman- 
Penrose formalism, for further details see Ref. [42]. That is, in the Kerr background we adopt the following vectors 



2 Even though the radial and angular equations are coupled via the separation constant, qAi m , we are able to find excellent agreement 
with the CFM by starting from the Schwarzschild (a = 0) result for our initial guess of the Kerr (a ^ 0) QNM solution using FindRoot 
in Mathematica in our AIM code (at least for a < 1). 
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TABLE VIII. The QNM Frequencies for the Scalar Perturbations of the Kerr Black Hole, with a = 0.80 (M = 1, m = 0). 
Numerical data via the CFM taken from [41], where the AIM was set to run at 15 iterations. 



1 
1 


n 


Numerical 


TV^Y-rl Oc^or WtTR/' "H 


oixtn L^rQer vvivjj^j j 




U 


U 


n 1 1 a c n nnK^ 

u.ii4o -u.uyo (i 


0.1005 - 0.100 II 


U.1211 - U.Uoy (1 


nil/11 n noon; 








/ 10 00/ K 000/ \ 


(0. /D/0 ,- O.ZO/oJ 


^-U.oo/o , -l.oo/oj 


1 





0.3067 -0.0901i 


0.3029 - 0.0891i 


0.3053 - 0.0893i 


0.3052 - 0.0892i 








(-1.2470 , -1.11%; 


(-0.4o7o , -0.89/0 ) 


/ n /inoy 1 nn07\ 
(-0.4970 , -I.OO70) 




1 


0.2820 -.2.1061 


C\ O l "7tT O n 077H " 

0.2758 - 0.2779i 


c\ 0001 c\ 07c 

0.2821 - 0.27501 


0.2817 - 0.2756i 










fn n4% 1 ni %1 




2 





0.5071 -0.0897i 


0.5035 - 0.0885i 


0.5041 - 0.0886i 


0.5041 - 0.0886i 








(-0.71%, -1.34%) 


(-0.59%, -1.23%) 


(-0.59% , -1.23%) 




1 


0.4906 -0.2722i 


0.4866 - 0.2693i 


0.4885 - 0.2690i 


0.4885 - 0.2689i 








(-0.82%, -1.07%) 


(-0.43%, -1.18%) 


(-0.43% , -1.21%) 




2 


0.4609 -0.4634i 


0.4585 - 0.4570i 


0.4607 - 0.4581i 


0.4606 - 0.4579i 








(-0.52%, -1.38%) 


(-0.04%, -1.14%) 


(-0.06%, -1.19%) 



as the null tetrad: 



lj = ^-(A,-p 2 ,0,-aAsin 2 i 



1 

-(«asin6>, 0, -p 2 , -i(r 2 + a 2 ) sin6>) , 



-(A,p 2 ,0,-aA sin 2 ' 

1 



J V2p 
V = l(r 2 + a 2 ,A,0,a) , 

' (r 2 + a 2 ,-A,0,a) , 



V2p 

m? = — ^(iasin#, 0, 1, — i ) , 

y/2p sinlT 

where fhj and fh? are nothing but complex conjugates of rrij and m? respectively. 

It is clear that the basis vectors basically become derivative operators when these are applied as tangent vectors to 
the function e l ( UJt + m 't>) , Therefore we can write 

1 = D = V ,n = D* = ,m = S= -^_C\ ,m = 6* = , (5.16) 

where, 

^ ^ K n r-M 
T> n = d r + i—+ 2n , 

^+ n K n r-M 

V\ = d r -i-+ 2n^^- , 

£ n = dg + Q + n cot 6 , 
= d g - Q + n cot 6 , 

and X = (r 2 + a 2 )w + am with Q = sin 9 + to esc (9. 

The spin coefficients can be written as a combination of basis vectors in the Newman-Penrose formalism which are 
now expressed in terms of the elements of different components of the Kerr metric. So by combining these different 
components of basis vectors in a suitable manner we get the spin coefficients as 

k = a = \ = v = e = . (5.17) 
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V2p 



1 cot 6» 

p* 2y/2p 

18 A 



sin t 



, I" = 



r - M 



2p*p 2 
a = tt — /3* . 



7 = p + 



2p 2 



(5.18) 



Using the above definitions, and by choosing: fx — p*F%, g 2 = p(*2, fi 
pair of spinors) the Dirac equation reduces to 

1 



=£4/2 = 0, 



API / 2 - ^t/i = 0, 



i*2 and gi = Gi (where F± t 2 and G1.2 are a 



(5.19) 



1 



\/2 2 







AV\gi+V2Cig 2 = 0. 

2 2 

We separate the Dirac equation into radial and angular parts by choosing, 

AM) = R_ h (r)S_i(ff) , h{r,0) - iii(r)Si(fl) , 
ffi(r,e) = fl J (r)S_ i (fl) ,g 2 (r,9) = R_ i (r)S i (6) . 
Replacing these fj and gj (J — 1, 2) and using A as the separation constant, we get, 

C1S1 — —\S_i , 



c{s_- 



XSi 



(5.20) 



A5X> i?_i = AA5i?i , 

2 2 

A5X>^A5i?i = XR_i , 



(5.21) 



where 2 2 i?_ 1 is redefined as R_ 1 . 

2 2 

Eqs. (5.20) and (5.21) are the angular and radial Dirac equation respectively, in a coupled form with the separation 
constant A [38]. Decoupling Eq. (5.20) gives the eigenvalue/angular equation for spin half particles as 



CiC\ + A 2 

and Si satisfies the 'adjoint' equation (obtained by replacing 9 by n 
radial equation for spin half particles as 



AV\V Q - A 2 



S_i=0, (5.22) 
. 3 Decoupling Eq. (5.21) then gives the 

R_t=0, (5.23) 



and A2 Ri satisfies the complex-conjugate equation. Furthermore, unlike the case of a scalar particle, a spin half par- 
ticle is not capable of extracting energy from a rotating black hole, that is, there is no Penrose Process (superradiance) 
equivalent scenario [42]. 

Returning now to the AIM, recall that it shall work better on a compact domain, where we define a new variable 
y 2 = 1 — which ranges from at the event horizon (r = r+) to 1 at spatial infinity. It is then necessary to 
incorporate the boundary conditions, which expressed in the new compact domain is 



4>(y) 



-(i-z/ 2 ) 



— 77 — 1(7- 



(y 2 )*- a +(i-y 2 r r+ "e MT ^x(y), 



(5.24) 



3 Note that this angular equation is that given by Eq. (5.4) for s = 1/2 and hence using the method in Appendix A, we could solve this 
numerically. 
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where we have defined 



and ip satisfies the WKB(J)-like equation: 



tP = VAR 1/2 + R_ 1/2 (5.25) 



+ (w 2 - V)tp = (5.26) 



d 2 ^ 
dy 2 

with potential 

and K = K/uj = (r 2 + a 2 ) + am/oj (for more details see [42]). 

By making the change of coordinates and change of functions, our equation takes the form 

x(y) = A (t/) + s (y), (5.28) 

where as in Sub-Sec. V A we have 

A„ = -2^-^, (5.29) 

gay j ay 

and 

1 d 2 5 1 df 1 dg 1 , , , , 

*° = ~ 7d£ x * *T 7* (w ~ ■ (5 - 30) 

As presented in Table IX and Table X the QNM frequencies for the spin half perturbations of the Kerr black hole 
with the two "extreme" values of the angular momentum per unit mass, that is a — 0.00 and a = 0.80, m was set to 
0, while I was given values of 0, 1 and 2 and n varied accordingly. 

Included in Table IX are the numerically determined QNM frequencies published by Jing et al. [43]. Even though 
the WKB method has been used to calculate the Schwarzschild limit QNM frequencies before [24], the sixth order 
WKB values and AIM values are novel to this work and shall be explored more fully in Ref. [42]. The percentages 
bracketed under each QNM frequency, is the percentage difference between the calculated value and the numerical 
value published by Jing et al. [43]. 

As expected, since there are additional correction terms, the sixth order WKB QNM frequencies are closer to the 
numerical values than the third order WKB values. While the AIM does not prove as accurate in its calculation of 
the spin half QNM frequencies as it did with the scalar values (both for 15 iterations), none of the differences between 
the AIM values and the numerical values exceed 0.30%, except for when I — 2, n — 2 (better accuracy can be achieved 
by increasing the number of iterations). 

Similarly in Table X are the numerically determined QNM frequencies published by Jing et al. [43]. Both the 
third and sixth order WKB values along with the AIM values are novel to this work and shall also be explored more 
fully in Ref. [42] . The percentages bracketed under each QNM frequency, are the percentage differences between the 
calculated value and the numerical value published by Jing et al., at least for I = and I — 1. For 1 = 2, the AIM 
values are compared to the sixth order WKB values. As already noted, since there are additional correction terms, 
the sixth order WKB QNM frequencies are closer to the numerical values than the third order WKB values. Again 
the AIM does not appear to be as precise in calculating the QNM frequencies for spin half perturbations of the Kerr 
black hole as it was for the scalar perturbations (at least for 15 iterations). As we mentioned, additional tables and 
plots of these Kerr processes shall constitute a future work [42] . 
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TABLE IX. The QNM Frequencies for the Spin Half Perturbations of the Kerr black hole, with a = 0.00, that is, the 
Schwarzschild limit (M = 1, m — 0). Numerical data via the CFM taken from [43], where the AIM was set to run at 15 
iterations. 



1 


n 


Numerical 


Third Order WKB 


Sixth Order WKB 


AIM 








0.1830 - 0.0970i 


0.1765 - O.lOOli 


0.1827 - 0.0949i 


0.1830 - 0.0969i 








(-3.55% , 3.20%) 


(-0.16% , -2.16%) 


(<0.01% , -0.10%) 


1 





0.3800 - 0.0964i 


0.3786 - 0.0965i 


0.3801 - 0.0964i 


0.3800 - 0.0964i 








(-0.37% , 0.10%) 


(0.03% , <0.01%) 


(<0.01% , <0.01%) 




1 


0.3558 - 0.2975i 


0.3536 - 0.2987i 


0.3559 - 0.2973i 


0.3568 - 0.2976i 








(-0.62% , 0.40%) 


(0.03% , -0.07%) 


(0.28% , 0.03%) 


2 





0.5741 - 0.0963i 


0.5737 - 0.0963i 


0.5741 - 0.0963i 


0.5741 - 0.0963i 








(-0.07% , <0.01%) 


(<0.01% , <0.01%) 


(<0.01% , <0.01%) 




1 


0.5570 - 0.2927i 


0.5562 - 0.2930i 


0.5570 - 0.2927i 


0.5573 - 0.2928i 








(-0.14% , 0.10%) 


(<0.01% , <0.01%) 


(0.05% , 0.03%) 




2 


0.5266 - 0.4997i 


0.5273 - 0.4972i 


0.5265 - 0.4997i 


0.5189 - 0.5213i 








(0.13% , -0.50%) 


(-0.02% , <0.01%) 


(-1.46% , 4.32%) 



TABLE X. The QNM Frequencies for the Spin Half Perturbations of the Kerr black hole, with a = 0.80 (M — 1, m = 0). 
Numerical data via the CFM taken from [43], where the AIM was set to run at 15 iterations. 



1 


n 


Numerical 


Third Order WKB(J) 


Sixth Order WKB(J) 


AIM 








0.1932 - 0.089H 


0.1883 - 0.0896i 


0.1914 - 0.0865i 


0.1920 - 0.0872i 








(-2.54% , 0.56%) 


(-0.93% , -2.92%) 


(-0.62% , -2.13%) 


1 





0.3993 - 0.0893i 


0.3956 - 0.088L 


0.3967 - 0.0880i 


0.3965 - 0.0880i 








(-0.93% , -1.34%) 


(-0.65% , -1.46%) 


(-0.70% , -1.46%) 




1 


0.3789 - 0.2728i 


0.3751 - 0.270L 


0.3777 - 0.2687i 


0.3764 - 0.2517i 








(-1.00% , -0.99%) 


(-0.32% , -1.50%) 


(-0.66%, -7.73%) 


2 







0.5984 - 0.088L 


0.5987 - 0.088H 


0.5987 - 0.0882i 












(<0.01% , 0.11%) 




1 




0.5844 - 0.2669i 


0.5855 - 0.2667i 


0.5846 - 0.2644i 












(-0.15% , -0.86%) 




2 




0.5600 - 0.4512i 


0.5609 - 0.4517i 


0.6023 - 0.4260i 












(7.38%, -5.70%) 



C. The Spin- Two Case 

As we have mentioned earlier, the radial equation for nonzero spin s is in general complex. In fact, it does not even 
reduce to the Regge- Wheeler and Zerilli equations when the rotation parameter a — > 0. Detweiler [40] has found a 
way to overcome this problem, where he defined a new function 



X 



A*/ 2 (r 2 +a 2 ) 1/2 



a{r)R + /3(r)A 



3+1 dR 
dr 



(5.31) 



If the functions a(r) and /3(r) are required to satisfy 



a 2 - a'[3A s+1 + a(3'A s+1 - fi 2 A 2s+1 K = constant , 



(5.32) 
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then it can be shown that the radial equation in Eq. (5.5) becomes 

d 2 X 



dx 2 



VX = 



where 



V 
G 



AU dG 
~ 2+G + dx~> 



(r 2 + a 2 ) 
s(2r- 1 



rA 



U = K- 
x = r + 



2{r 2 + a 2 ) ( r 2 + a 2 ) 2 ' 
2a' + (P'A S+1 )' 



PA S 

— — In (r — r + ) In (r — r_) 

r _ r _ r _ r _ 



(5.33) 

(5.34) 
(5.35) 

(5.36) 
(5.37) 



As Detweiler has indicated, it is possible to choose the functions a(r) and /3(r) so that the resulting effective potential 
V(r) is real and has the form 

P 2 A ' 



V 



{r 2 + a 2 f 1 9 + bA 



/(/ + 2) bA | [KpA-(g'A-gA')][ K pg-b(g'A-gA')] 



p 2 (g + bA) (g-bAY 



ramA 



top (r 2 + a 2 Y 



d 



f 2 + a 2 dr 



ramA 



uip (r 2 + a 2 ) 



Y 2 



(r 2 + a 2 Y 



where 



g = fp 2 + + a 2 ) - ir 2 A , 

9 , am 
p = r 2 + a 2 , 



K = ± 



{9 - 2/ [(a 2 ~ ^) (5/ + 6) - 12a 2 ] + 2bf(f + 2)} 



1/2 



a 

am\ 



^ ± 3(« 2 -- y , 
Y = am — (r 2 + a 2 )u> , 
f = A + a 2 uj 2 — 2amoj 



(5.38) 

(5.39) 
(5.40) 

(5.41) 

(5.42) 
(5.43) 
(5.44) 



When the Kerr rotation parameter a approaches zero, the potential V in Eq. (5.38) coincides with the Regge- Wheeler 
potential for negative k, and coincides with the Zerilli potential for positive n. Here we choose k to be negative, where 
the choice of the sign in Eq. (5.42) is determined by the sign of m [7, 40]. 
The QNM boundary conditions for X are 



X 



e ikx 



(5.45) 



where 



Hence, we write 



k = lj 



X = e iujr r r +~ r - 



(r-r + ) r + 
(r — r_) r - 



Substituting this into Eq. (5.33), we have 



XG,: 



2IY 



A(r 2 + a 2 ) 



XG,r + 



Xg 



+ ^G,r 



T G (r 2 - a 2 
A(r 2 +a 2 ) 



2 1 ? \ ^ 
r* + a z 



A 



V 



Xg = , 



(5.46) 



(5.47) 



(5.48) 



2G 



TABLE XI. Spin-2 angular separation constants and Kerr gravitational quasinormal frequencies for the fundamental mode 
corresponding to I = 2 and m = compared with the CFM [6] (M = 1/2). 



a 


-^Leaver 


Aaiai 


(-d Leaver 







(4.0000, 0.0000) 


(4.0000, 0.0000) 
(<0.01%)(<0.01%) 


(0.7473, -0.1779) 


(0.7413, -0.1780) 
(-0.80%) (-0.06%) 


0.1 


(3.9972, 0.0014) 


(3.9973, 0.0014) 
(<0.01%)(<0.01%) 


(0.7502, -0.1774) 


(0.7444, -0.1775) 
(-0.77% , -0.06%) 


0.2 


(3.9886, 0.0056) 


(3.9887, 0.0056) 
(<0.01%)(<0.01%) 


(0.7594, -0.1757) 


(0.7540, -0.1763) 
(-0.71%) (-0.03%) 


0.3 


(3.9730, 0.0126) 


(3.9733, 0.0126) 
(<0.01%)(<0.01%) 


(0.7761, -0.1720) 


(0.7715, -0.1722) 
(-0.59%) (-0.12%) 


0.4 


(3.9480, 0.0223) 


(3.9482, 0.0222) 
(<0.01%)(-0.45%) 


(0.8038, -0.1643) 


(0.8025, -0.1639) 
(-0.16%) (0.24%) 


0.45 


(3.9304, 0.0276) 


(3.9303, 0.0280) 
(<0.01%)(1.45%) 


(0.8240, -0.1570) 


(0.8250, -0.1591) 
(0.12%)(-1.34%) 



where 



rv 



iw + 



ICJ 



r(r + — r_) 



ikr 



(5.49) 



As we did earlier, we define the variable £ = 1 — r + /r which has a compact domain < £ < 1. Eq. (5.48) can then 
be written in the AIM form 



Xgm = Ag(£)xg,« + s g{£)xg 



where 



i-£ (i-£) 2 



2r 



1 H 



G 



a 2 (l-0 S 



A 



•a 2 (l-0 s 



r r 2 +fl 2 (1 _ 2 
A(l-0 2 



(i-0 4 



r G ^ - a 2 (i-fl- 

A r5.+o a (l-f) s 



(i-O s 



(5.50) 

(5.51) 
(5.52) 



The results for the gravitational (spin-two) case are presented in Tables XI and XII. In general the error in the 
separation constant is smaller than that of the quasinormal frequencies. As for the quasinormal frequencies the error 
in the Kerr case is larger than that of either the Schwarzschild or the Reissner-Nordstrom cases, where this is due 
to our consideration of the angular and the radial equations simultaneously. The number of iterations that can be 
performed in the code is relatively small, much like the number of continued fractions in the CFM is typically smaller 
due to the coupling between radial and angular equations. 



VI. DOUBLY ROTATING KERR (A)DS BLACK HOLES 

Rotating black holes in higher dimensions were first discussed in the seminal paper by Myers and Perry [44] . One of 
the unexpected results to come from this work was that some families of solutions were shown to have event horizons 
for arbitrarily large values of their rotation parameters. The stability of such black holes is certainly in question 
[45, 47], with numerical evidence recently provided by Shibata and Yoshino [46]. 

Another new feature of the Myers-Perry (MP) solutions is that they in general have L^i~^J spin parameters, making 
them more complex than the four dimensional Kerr solution. The first asymptotically non-flat five-dimensional MP 
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TABLE XII. Spin-2 angular separation constants and Kerr gravitational quasinormal frequencies for the fundamental mode 
corresponding to I = 2 and m = 1 compared with the CFM [6] (M — 1/2). 



a 


-^-Leaver 


A AIM 


^Leaver 


maim 





(4.0000, 0.0000) 


(4.0000, 0.0000) 
(<0.01%)(<0.01%) 


(0.7473, -0.1779) 


(0.7413, -0.1780) 
(-0.80%) (-0.06%) 


0.1 


(3.8932, 0.0252) 


(3.8937, 0.0250) 
(-0.01%) (-0.89%) 


(0.7765, -0.1770) 


(0.7726, -0.1755) 
(-0.51% , 0.86%) 


0.2 


(3.7676, 0.0532) 


(3.7681, 0.0526) 
(0.01%)(-1.12%) 


(0.8160, -0.1745) 


(0.8143, -0.1726) 
(-0.02%) (1.09%) 


0.3 


(3.6125, 0.0835) 


(3.6123, 0.0826) 
(<0.01%)(-0.99%) 


(0.8719, -0.1693) 


(0.8722, -0.1674) 
(0.03%)(1.00%) 


0.4 


(3.4023, 0.1122) 


(3.4011, 0.1110) 
(-0.03%) (-1.00%) 


(0.9605, -0.1559) 


(0.9620, -0.1543) 
(0.16%)(1.05%) 


0.45 


(3.2535, 0.1195) 


(3.2491, 0.1173) 
(-0.03%) (-1.82%) 


(1.0326, -0.1396) 


(1.0376, -0.1369) 
(0.48%) (1.97%) 



metric was given in Ref. [48]. Subsequent generalizations to arbitrary dimensions was done in Ref. [49], and finally 
the most general Kerr-(A)dS-NUT metric was found by Chen, Lii and Pope [50]. 

In this section we review how the AIM can be used to solve the D > 6 two-rotation scalar perturbation equations 
(for more details on the metric and resulting separation see Ref. [51]). The scalar field master equations are found to 
be [51]: 
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(6.1) 

(6.2) 

(6.3) 
(6.4) 

(6.5) 

(6.6) 
, and 



and i = 1,2. In the above g is the curvature of the spacetime satisfying R^ = —'ig 2 g^ v (e.g., see Ref. [50] 
ax, a,2 are the two rotation parameters and for later reference we define e = aifa\. 

Doubly rotating black holes are more complicated than simply rotating black holes (cf. Ref. [52]), because two 
rotation planes lead to two coupled spheroids which are also needed for the solution of the radial equation. 



A. Radial Quasi-Normal Modes 



For simplicity we will consider the flat case, setting g = 0, which leads to easier QNM boundary conditions (cf. 
Schwarzschild to Schwarzschild-dS). These satisfy the boundary condition that there are only waves ingoing at the 
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black hole horizon and outgoing waves at asymptotic infinity. 

As we have shown with the previous examples, it is easier to work on a compact domain and define the variable 
x = 1/r, so that infinity is mapped to zero and the outer horizon stays at Xh — = 1. The domain of x will 

therefore be [0, 1]. Thus the QNM boundary condition is translated into the statement that the waves move leftward 
at x = and rightward at x = 1. We again choose the AIM point in the middle of the domain, that is, at x — 1/2. 

In terms of x the radial equation (6.1) becomes: 

= -x°~*l (-*-"A.g) + ( ( ^ 2 + a ^f 2 + ai)2 ^ - alaUU + D - 7)x* - | - 6 2 ) R , (6.7) 

where A x (a;) = A r (r = l/x) and u> x (x) = co r (r = 1/x). 

After performing some asymptotic analysis, we find that for the solutions to satisfy the QNM boundary conditions 
we must have: 

R ~ (1 - x)^ hah x {D - 2)/2 e l ^ /x y{x) , (6.8) 

where 

LU h = U] X (X = 1) , 

(l + qj)(l + al) 
" _ Ai(x = l) ' 
We then substitute this ansatz into Eq. (6.7) and rewrite into the AIM form 

y" = X y' + s Q y . (6.11) 

This final step above can be performed in Mathematica, where the resulting expressions for Ao and sq are fed into the 
AIM routine. The method we use to find the QNMs proceeds in a fashion similar to that used in Ref. [53, 54] (see 
also Sec V C) except we use the AIM instead of the CFM. 

First we set the number of AIM iterations in both the eigenvalue and QNM calculations to sixteen. We start with 
the Schwarzschild values (01,02, w), that is, at the point (01,02) ~ and then increment ai and 02 by some small 
value. We take the initial eigenvalues (61, 62), insert them into the radial equation (6.8) then use the AIM to find the 
new QNM that is closest to u using the Mathematica routine FindRoot. 

Taking this new value of omega, u/, we insert it into the two angular equations (at the same value of a\ and 02) 
then solve using the AIM, searching closest to the previous b\ and 62 values. Thereby obtaining the new eigenvalues 
b'i, b' 2 - We then repeat this process with the new (u/,&i,b 2 ) as the starting point until the results converge and we 
have achieved four decimal places of accuracy. When this occurs we increment a\ and 02 again and repeat the process. 
In this way, we are able to find the QNMs and eigenvalues along lines passing approximately through the origin (that 
is, starting from the near Scwharzschild values) in the (01,02) parameter space. 

As an example, we have plotted various values of e = 02/01 (= 0,0.2,0.4,0.6,0.8,1) against a\ and used an 
interpolating function to interpolate between these values as shown in Fig. 1 (for further details see Ref. [51]). 

VII. SUMMARY AND OUTLOOK 

In this review we have shown that the AIM can be used to calculate the radial QNMs of a variety of black hole 
spacetimes. In particular, we have used it to calculate perturbations of Schwarzschild (in asymptotically flat, de Sitter 
and anti-DeSitter), RN and Kerr (for spin 0, 1 and 2 perturbations) black holes in four dimensions. We argued that 
the method will be of use in studies of extra dimensional black holes and gave an explicit example of this in the case 
of the doubly rotating Myers Perry black hole. 



(6.9) 
(6.10) 
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FIG. 1. An example of the D = 6 fundamental (J, mi, m2, rii, n.2) = (0,0,0,0,0) QNM. On the left is a plot of the imaginary 
part and on the right plot of the real part. 

We have hopefully demonstrated how the AIM can also be applied to radial QNMs and not just to spheroidal 
eigenvalue problems [11, 31]. Given the fact that the AIM can be used in both the radial and angular wave equations 
[31] we expect no problems in obtaining QNMs for Kerr-dS black holes in four and higher dimensions. Note that 
this was only recently accomplished via the CFM in Ref. [55] using Heun's equation [56] to reduce the problem to a 
3-term recurrence relation. In higher dimensions a similar method was used for simply rotating Kerr-AdS black holes 



It remains to be seen if the AIM can be tailored to handle asymptotic QNMs (see Ref. [57] for an adapted version 
of the CFM) . However, given the close relationship between the AIM and the exact WKB approach [58] , it might be 
possible to adapt the AIM to find asymptotic QNMs [59-62] numerically or even semi-analytically. 

The AIM might be of some topical use, for example, in the angular spheroids/QNMs needed in the phenomenology 
of Hawking radiation from spinning highcr-dimcnsional black holes, for a recent review see Ref. [63]. We recently used 
a combination of all the techniques discussed in this work to evaluate the angular eigenvalues, o^kjmi f° r rea l c = aw, 
which are needed for the tensor graviton emission rates on a simply rotating Kerr-de Sitter black hole background 
in (n + 4)-dimensions [64] (also see Ref. [65]) and it might also be interesting to find QNMs of doubly rotating 
Kerr-(A)dS black holes (for asymptotically flat Kerr see Ref. [51]). Finally, attempting to solve the QNMs for all 
spins on the Schwarzschild-AdS background via the AIM also seems an interesting problem. 

As such we hope to have provided the reader with enough technical details, and to have addressed some of the 
possible questions to allow them to pursue the study of QNMs with the AIM. 4 
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4 Source code and other information for some of the cases presented here can be found on the AIM link at http : //www-het . phys . sci . 
osaka-u. ac . jp/~naylor/. 
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TABLE XIII. Selected spin two eigenvalues, iAi m , obtained from the AIM for a Kerr black hole with varying values of c= aui 
and m = 1. The same number of iterations in the AIM, n A and the number of recursions in the CFM, nc for c — 0.1, 0.8 at a 
working precision of 15 digit precision, where results are presented to 10 s.f. 



1 


c = 0.1 (n A =n c = 15, ) 


c = 0.8 (n A =35, n c =70) 


c = -lOi (n A =80, n c =130 CFM) 


c = 10 (n A =75, n c =145) 


2 


-0.1391483511 


-1.462479552 


(12.44128209, 0.8956143162) 


-101.8949078 


3 


5.929826236 


5.247141863 


(32.31138608, 1.302608040) 


-63.74900642 


4 


13.95640426 


13.45636668 


(51.27922784, 1.946041848) 


-30.35607486 


5 


23.96944247 


23.54163307 


(69.25750923, 3.012877154) 


-4.557015739 


6 


35.97681567 


35.58524928 


(85.86796852, 4.990079008) 


-2.555206382 


7 


9.98139515 


49.61077286 


(99.20081385, 6.801617108) 


12.32203552 
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Appendix A: Angular Eigenvalues for Spin- Weighted Spheroidal Harmonics 



As mentioned in Sec. VC, aside from radial QNMs the AIM can also be applied to various kinds of spin- weighted 
spheroidal harmonics, s Si m (0), e.g. see Ref. [54]. Therefore, in this appendix, we briefly compare the AIM with the 
CFM for the four-dimensional spin-weighted spheroids. 

With the regular boundary conditions, the angular wave function can be written as [6] 



S = e au)u (l + u)^ m - sl (l - u)^ m+s \ XA ( u ) . 
Putting this back into Eq. (5.4) and rewriting the equation in AIM form, we have 

XA,uu = Xa(u)Xa,u + sa{u)xa , 



where 



\ A (u) 
s A (u) 



2u 



l-u 2 

1 



l-u 2 
N = ulu + 



- 2N 
(to + 



l-u 2 
\m — s\ 
2(1 + u) ~ 



+ 2au)su + 2uN — (a 2 uj 2 u 2 + s + s Ai m ) 



N — N , 



\m ■ 



u) 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 



These are the relevant equations for calculating the eigenvalues of the spheroidal harmonics in the four-dimensional 
case. It was noticed in Ref. [11] that the AIM converges fastest at the maximum of the potential, when the AIM is 
written in WKB form. In four dimensions this occurs at x — and is true for general spin-s as we have verified. Note 
that for higher dimensional generalizations it is not easy to explicitly find a maximum [31]. 

It may be worth mentioning that for the case where c = an exact analytic solution of the above equations leads 
to [54]: 



,Ai m = 1(1 + 1) - s(s + 1) 



(A6) 



In the exact limit c = the AIM recovers the result for spherical harmonics, Eq. (A6) above, while for the CFM 
taking c = leads to singularities [54]; however for c <C 1 we find agreement with the CFM and Eq. (A6). 

For the purposes of consistency we have calculated (see Table XIII) the 2^4;i eigenvalues for the lowest I — 2, ... ,7 
modes to 10 significant figures and have also compared this with the results of the CFM. In both the AIM and 
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FIG. 2. (Color Online) D — 6, g = 0, (j, mi, mj, ru, ni) = (0,1,1,0,0). A plot of the eigenvalues for various choices of 
e = aija\. Note that the dependence on a\ has been scaled into the other quantities. 




FIG. 3. (Color Online) D — 6, e = 0,2/0,1 = 1/2, (j, mi, TO2, Tii, 712) = (0, 1, 1, 0, 0). A plot of the eigenvalues for ga\ = 0.5i, 0, 0.5, 
corresponding to deSitter, flat, and anti-deSitter spacetimes respectively. Note that the dependence on a\ has been scaled into 
the other quantities. 

CFM larger I modes require more iterations/recursions to achieve convergence in a given I eigenvalue to the required 
precision. Care should be taken when comparing the number of iterations in the AIM with that of the number of 
recursions in the CFM, because one iteration of the AIM is not equivalent to one iteration of the CFM. In fact 
although we typically need to iterate the improved AIM on average a lesser number of times, the CFM is typically 
faster for smaller values of c. However, for larger values of c both methods can be faster or slower. 

The results of the first few I eigenvalues for different values of c = au, with m = 1, are presented in Table XIII. As 
far as we are aware this is the first time tables of spin-2 spheroids (for general complex c) have been presented using 
the AIM. Further results are presented in Sec. V C along with the radial QNMs for the spin two perturbations of the 
Kerr black hole in Table XII. 



Appendix B: Higher Dimensional Scalar Spheroidal Harmonics with two Rotation Parameters 

The two Eqs. (6.2) are in fact the two-rotation generalization of the higher dimensional spheroidal harmonics 
(HSHs) studied in Ref. [53]. In this case, the existence of two rotation parameters leads to a system of two coupled 



32 



second order ODEs 5 . In general, one would expect that the generalizations of the HSHs to L^-^J rotation parameters 
would lead to even larger systems of equations. While these systems would also be useful generally in studies of MP 
black holes, here we will only focus on the two rotation case. 

The angular equations can be written in the Sturm-Liouville form (assuming momentarily that ui and 62 are real): 

Atu = -^r (?>(&) (&)) +Q^i)ReM (Bl) 

with the weight function wi(£i) = 5 ^ 2 , the eigenvalue A = —b\, and 

P«i) = «, <£ "" 2 4 { ,. (B2) 

m _ i ?f - w , ^-(.m-u' q + _ ^ (B3) 

where and uj^ are defined in the obvious way under the change of coordinates. Since w(£) > we can define the 
two norm's: 

^(iZjJa r ^ D - 5)/2 \R 01 \ 2 dCi , (B4) 

a 2 

N 2 (Re 2 )* [ \i D - 5)/2 \R e . 2 \ 2 d^ . (B5) 
Jo 

For further details see Ref. [51]. 

The regular solutions are found to be: 

JZi~ (£1-0^(0? -6)^*1; &€(al,a?), (B6) 
fl2~3 /2 (ai- 6)^*2; &e(o,a£). (B7) 

Now for a given value of u we can determine b\ and 62 simply by performing the improved AIM [10] on both of 
the angular equations separately. This will result in two equations in the two unknowns b\ , hi which we can then be 
solved using a numerical routine such as the built-in Mathematica functions NSolve or FindRoot. More specifically 
we rewrite Eqs. (6.2) using (B6) and (B7) and transform them into the AIM form: 

-w = Ao1 k + -*■• (B8) 

-T7T = A 02 — - + S 2*2 • B9 
d^i d& 

The AIM requires that a special point be taken about which the Aoi and soi coefficients are expanded. As was shown 
in Ref. [31] different choices of this point can worsen or improve the speed of the convergence. In the absence of a 
clear selection criterion we simply choose this point conveniently in the middle of the domains: 

&i = ^. te-f (BIO) 

Some results are plotted in Figs. 2, 3 above. This method can now be used in the radial QNM equation in Sec. 
VIA. 
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